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ABSTRACT 

Aims. We investigate the long-term dynamical evolution of two distinct stellar populations of low-mass stars in globular clusters in 
order to study whether the energy equipartition process can explain the high number of stars harbouring abundance anomalies seen in 
globular clusters. 

Methods. We analyse N-body models by artificially dividing the low-mass stars (in < 0.9 Mq) into two populations: a small number 
of stars (second generation) consistent with an invariant IMF and with low specific energies initially concentrated towards the cluster- 
centre mimic stars with abundance anomalies. These stars form from the slow winds of fast-rotating massive stars. The main part of 
low-mass (first generation) stars has the pristine composition of the cluster We study in detail how the two populations evolve under 
the influence of two-body relaxation and the tidal forces due to the host galaxy. 

Results. Stars with low specific energy initially concentrated toward the cluster centre need about two relaxation times to achieve 
a complete homogenisation throughout the cluster For realistic globular clusters, the number ratio between the two populations 
increases only by a factor 2.5 due to the preferential evaporation of the population of outlying first generation stars. We also find that 
the loss of information on the stellar orbital angular momentum occurs on the same timescaie as spatial homogenisation. 
Conclusions. To reproduce the high number of chemically anomalous stars in globular clusters by preserving an invariant IMF, more 
efficient mechanisms such as primordial gas expulsion are needed to expel the stars in the outer cluster parts on a short timescaie. 
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' 1. Introduction 

In a given globular cluster (GC), low-mass stars are known to 
display a homogeneous composition in Fe and other heavy el- 
ements (with the noticeable exception of a> Cen). However, in 
each individual GC some pronounced star-to-star variations in 
light element abundances (Li to Al) are found, leading to the 
well-known C-N, O-Na, Li- Na and Mg-Al anticoiTelations (for 
a review see Gratton et alj2004; Charbonnel 2005) . They are the 
I direct results of H-burning nr i cleosynthesis at high temperatures 
' ('Denisenkov & Denisenkova] Il989l 1 19901: iLanger & HoffmanI 
il995; Prantzos et al. 2007). In addition, since such abundance 
patterns are not detected in population II halo fields stars, they 
are expected to be related to the evolution of the GCs. Since the 
low-mass stars still on the main sequence or on the red giant 
branch display such abundance variations and since these stars 
do not produce in their interior the high temperatures needed 
to create the abundance anomalies, it is expected that low-mass 
stars have inherited this chemical pattern at birth. 

Massive stars as possible sources to contaminate the 
metal abundances in the GCs have first been proposed by 
iPilachowski et al. ( 19821 ) and [W allerstein et al. (1 983) a r id the n 
studied by Prantzos & Charbonnel| (|2006|) and ISrnithl (l2006h . 
Recentlv Decressin et aL (2007bi) explore the role of a high ini- 
tial rotation in such stars, finding that abundance anomalies in 
GCs can result from an early pollution by fast-rotating massive 
stars. Indeed during main sequence evolution, angular momen- 
tum is transported from the centre to the stellar surface, and for 
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Stars heavier than 20 Mq with a high initial rotational velocity, 
their surface can reach the break-up velocity limit at the equa- 
tor (i.e. the centrifugal equatorial force balances gravity). In 
such a situation, a slow mechanical wind develops at the equator 
and f orms a disc around the stars similar to what h appens to Be 
stars dTownsend et all 120041: [Ekstrom et alj|2008l) . The second 
effect of rotation is to transport elements from the convective H- 
burning core to the stellar surface, which enriches the disc with 
H-burning material. While matter in discs has a slow outward 
velocity and hence will stay in the potential well of the cluster, 
matter released during the main part of the He-burning phase and 
during SN explosions has a very high radial velocity and will be 
lost from the cluster. Therefore, new stars can only form from 
the matter available in discs and can become the stars with abun- 
dance anomalies we observe today. Thus GCs can contain two 
distinct populations of low-mass stars: a first generation that has 
the chemical composition of the material out of which the clus- 
ter formed (similar to field stars with a similar metallicity); and a 
second generation of stars harbouring the abundance anomalies 
born from the ejecta of fast-rotating massive stars. 

Based on the determination of the c omposition of giant stars 
in NGC 6752 by I Carretta et all (l2007l) . iDecressin et al] ( l2007ah 
determined that around 85% of stars (of the sample of 120 stars) 
display abundance anomalies. Prantzos & Charbonnell (l2006h 
fi nd similar resul t s for N GC 2808 with their analysis of the data 
of ICarretta et al.l (l2006l) : 70% of the cluster stars present abun- 
dance anomalies. 

A main problem is that by assuming a lSalpeterl (Il955h IMF, 
the accumulated mass of the slow winds ejected by the fast- 
rotating massive stars would only represent 2.5% of the mass 
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of the whole populatio n of first generation stars (see Fig. 5 in 
iDecressin et al.ll2()07ab . This number is only a lower limit since 
some mixing between the slow wind and the gas of the interstel- 
lar medium left after the formation of first generation stars needs 
to be taken into account so that more matter will be available to 
form sec ond generation s tars. Ta king this dilution process into 
account, 'Decressi n et al.l (l2007ah found that the initial mass of 
the second generation stars is 5% of the cluster mass. As the 
whole population of low-mass stars represents half of the cluster 
stellar mass, this leads to a number ratio of about 0. 1 between the 
second and the first generation of low-mass stars. This dilution 
process of the slow wind in the ISM is also mandatory to ex- 
plain the Li det ected in stars w i th abu ndance anomalies. Indeed 
in NGC 6752, iPasquini et"!!] ( l2005h found that turn-off stars 
with abundance anomalies display a noticeable abundan ce vari- 
ation in Li. Similar results are found for 47 Tuc (Bonifaci o et alj 
I2007h . However in the hot interior of massive stars, all the Li 
initially present is consumed and no Li is ejected into the disc. 
Thus the Li in stars with abundance anomalies can only come 
from the ISM. 

To reconcile the high fraction of secon d generation stars with 
the observations IDecressin et al.l (l2007al) suggest that we need 
either (a) a flat IMF for the massive stars, or (b) that 96% of the 
low-mass stars of the first generation have escaped the cluster. 
The second case could be achieved if the following scheme hap- 
pens: we assume that the fast-rotating massive stars are born near 
the cluster centre or migrate rapidly toward the centre (i.e. , that 
the cluster is initially mass segregated). In such a case, the slow 
winds ejected by these massive stars are also concentrated near 
the centre. As the second generation stars need to be formed on 
a short timescale to avoid the ensuring supernovae pushing the 
slow wind material, from which second generation stars are to 
form, out of the cluster, we can expect that second generation 
stars will form with the same radial distribution as the massive 
stars. 

Primordial mass segregation is observed in young star 
clusters in the Milky Way and the Magellanic Clouds 
jBonnell & DaviesI 119981: iGouliermis eTaTI l2004t IChen et al.l 
I2007h . Thus we could reasonably expect that GCs were also 
mass segregated at birth. It should be noted that the recent corre- 
lation between the slope of the mass function and the concentra- 
tion parameter (c = log(r,/rc), where r t and rc are the tidal an d 
core radius respectively) discovered bv lDe Marchi et af] (l2007l) . 
can be explained n aturall y if the GCs star ted mass segregated 
jMarks et al.ll2008l:lBaumgardt et al.ll2008al) . 

Empirical evidence for the presence of gas in the cluster cen- 
tre at the same tim e as m assive stars are still evolving is provided 
by iGalliano et al.l (12008 ). They investigated extremely massive 
young clusters, which are believed to be nascent GCs (around 
6-8 Myr old), in the spiral Galaxy NGC 1365, and found an op- 
tically thick component in addition to a thin one. These authors 
argue that this thick component could be related to a subsequent 
or on-going episode of star formation. 

We expect that the slow winds of fast-rotating massive stars 
dilute the ISM locally around each massive star, so second gen- 
eration stars will be born only in the central part of the cluster 
while the outer regions are populated exclusively by first gener- 
ation stars. If some dynamical process unbinds the external parts 
most of the first generation stars may be lost while the second 
generation stars are retained in the cluster. In the present paper 
we investigate if the tidal interactions with the host Galaxy can 
produce such an effect. 

The time needed for the second generation of stars to spread 
throughout the cluster is expected to scale with the two-body 



rel axation time a t the half-mass radius, t±, of the cluster given 
bv lSpitzeii([T987l) . 
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where A^, and (m) are the number of stars and their mean mass 
while Mci, and rh represent the total mass and the half-mass ra- 
dius of the cluster respectively. The two-body relaxation time is 
an estimate of the timescale over which the cluster rearranges 
itself significantly due to its evolution towards energy equiparti- 
tion as driven by the many encounters between two stars within 
the cluster potential. At the same time, the tidal forces of the host 
galaxy remove stars from the outer cluster parts. Current obser- 
vations of GCs suggest that most of them are several relaxation 
times old, so that we would expect that stars in these clusters 
are thoroughly mixed and therefore homogeneous at the present 
time. 

In this paper, we study more quantitatively the interplay be- 
tween the effect of internal dynamics and the external potential 
of the host galaxy on the evolution of two stellar populations in 
GCs. In § 2, we describe the dynamical models we used. Then 
our analysis is presented in § 3. Finally, our results are discussed 
and our conclusions are drawn in § 4. 



2. Dynamical models 

2.1. Used models 

The analysis presented here is based on the calculations of 
Bau mgardt & Makinol (l2003h . which consist of a set of mod- 
els containing 8k to 128k star^ initially, c arried out wit h 
the collisional Aarseth N-body code nbody4 dAarsethl Il999l) . 
It uses a Hermite scheme with individual time-steps for 
the integration and treats cl ose en counters between stars by 
KS (Kustaanheimo & Stief el | 1965h and chain regularization 
(iMTkkola & Aarseth. 199a 1 1993^ 

Here we restrain the analysis to clusters moving on circu- 
lar orbits (with orbits 8500 pc away from the Galactic centre) 
through a host galaxy modelled as a logarithmic potential. The 
clusters follow King profiles initially with central concentrations 
W() - 5.0. The cluster initial radii are adjusted such that the 
tidal radius of the King model was equal to the tidal radius given 
by the host Galactic tidal field. The clusters are computed un- 
til complete dissolution. The initial mass function of stars in the 
clust er is given by a canonical two-part power-law mass func- 
tion (lKroupall200r ) with upper and lower mass limits of 15 
and 0.1 Mq respectively. This leads to an initial mean mass of 
(m) - 0.5 47 Mq. Stella r evolution is modelled by the fitting 
formulae of iHurlev et all (I2OOO1) with an assumed metallicity of 
Z = 0.001 ([Fe/H] ^ -1.2, flie mean value of GC metaUicities) 
for the cluster stars. Mass lost from the stars in the course of 
their evolution is assumed to be immediately lost by the cluster. 
The clusters contained no initial binaries. 



2.2. Selection by specific energy 

Although the models of lBaumgardt & Makinol (l2003l) have been 
computed for a single stellar population, we apply a process 
which could mimic the formatio n of a clus t er with two dynam- 
ically distinct populations. Decressin et al.l (l2007ah have shown 
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that the slow equatorial winds of fast-rotating massive stars cre- 
ate slowly expanding discs which stay in the potential well of 
the GCs. If we assume a radial velocity for these winds between 
1 and 10 kms"', and that the gas travels for about 10^ yr before 
starting to form new stars, it will travel for a distance of about 
0. 1 to 1 pc. As the initial tidal radius of the cluster is a few 10 pc, 
we expect that second-generation stars have the same initial spa- 
tial distribution as their massive progenitors. Since we assume 
here that the massive stars are born in the centre of the cluster, 
the second generation stars will also be located near its centre. 

In order to select stars in the cluster centre we sort all the 
low-mass stars (M < 0.9 Mq) according to their specific energy 
defined as the energy per unit mass. We define the second stellar 
generation as the stars with smallest specific energies, {i.e. those 
stars which are most tightly bound to the cluster due to their 
small central distance and low velocity). The number of second 
generation stars is given by having their total mass represent 5% 
of the initial mass of the cluster For the models with initially 
128k stars, we have 104691 first generation stars and 10327 sec- 
ond generation stars; the remaining stars are more massive than 

0.9 Mq. 

2.3. Possible issues 

It should be noted that the upper-limit for the IMF in the models 
o f Baumgardt & Makino (.2003) is 15 Mq. However the pollution 
by r nassive stars would invo lve stars more massive than 20 Mq 
(see iDecressin et al. I l2007bl) . The evolution of such high mass 
stars takes a few million years. So if the delay for the formation 
of the second generation stars is short, the two populations will 
be nearly coeval. However the direct dynamical effects of the 
fast-rotating massive stars (more massive than 20 Mq), like the 
lowering of the cluster potential well by the ejection of their fast 
polar winds and SN, or the creation of a second generation of 
stars with the matter released in their slow winds, are not taken 
into account in this paper. 

With the scheme presented in § 2.2, we implicitly assume 
that the lifetime of the first and second generation of low-mass 
stars are identical. In real clusters this would not be valid if 
a noticeable fraction of second generation stars shows a large 
overabundance in He. Indeed as He is the main product of H- 
burning in stars, abundance anomalies are also expected to come 
with He overabundance. Although He abundance cannot be di- 
rectly determined, several indirect observations push toward this 
finding. First, the GCs ai Cen and NGC 2 808 display multiple 
main sequences d Piotto et"al] I2005L l2007l): a doub le sub-giant 
branch is also found in NGC 1851 (iMilone etakl 12008). Such 
features can only be understood if the stars present various He 
contents. He enrichment is also a possible explanation for the 
appearance of extreme horizont al branches seen for several GCs 
(ICaloi & D'Antonal2005Ll2007h . 

For our paper, the different initial He content between the 
two populations needs to be properly taken into account as stars 
with higher initial He content have a shorter lifetime due to 
faster stellar evolution. Therefore to be consistent with the evo- 
lution of the second generation of stars requires custom-made 
models. However the value of F = 0.4 (instead of 0.25 ) is of- 
ten cited as needed to p roperl y explain the observation s jNorrisI 
I2004t iPiotto et al.ll2005h while IDecressin et all (l2007al) find that 
fast-rotating massive stars could even create low-mass stars with 
a He mass fr action up to F = 0.7. For 12 Gyr old clusters, 
iD'Antona & Ventura (2007) show that the turn-off mass de- 
creases from 0.8 to 0.64 Mq when the initial He content raises 
from 0.24 to 0.4. It could even be as low as 0.35 Mq if the ini- 
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Fig. 1. Position of low-mass stars (Mini < 0.9 Mq) in the Z- 
X (left) and Z-Y (right) plane at different times: Gyr (top), 
5.06 Gyr (middle) and 11.85 Gyr (bottom) for the N-body in- 
tegration with initially 128k stars. For sake of clarity, the upper 
and bottom parts of each panel display only the second and first 
generation of stars respectively. 



tial He content is 0.7 (Decre ssin et al., in prep.). Unfortun ately 
the N-body computations of iBaumgardt & Makinol (l2003h use 
only standard He abundance for the stellar evolution models, so 
this effect cannot be investigated in the present paper However 
we expect that the global findings of this paper on the dynam- 
ical timescales will not be changed by this He-rich population. 
Nevertheless a forthcoming paper will address this issue. 



3. Analysis 

3. 1 . Dynamical evolution: radial distribution 

Fig.[T]illustrates the evolution of two selected stellar populations 
in a cluster with initially 128k particles in a circular orbit around 
the Galaxy at three different times. In Fig. |2] one can see the 
radial distribution of the two populations at the same epochs. 
Initially, the second generation stars with low specific energy are 
concentrated within 6 pc around the centre of the cluster while 
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Fig. 2. Radial distribution of the first (dashed lines) and second 
(full lines) generation of low-mass stars at the same times as in 
Fig. [1] Each histogram is normalised to the total number of stars 
in each population. 
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Fig. 3. Evolution of the half-mass radius for first (dashed lines) 
and second (full lines) generation of low-mass stars present in 
the clusters as a function of time for the models with initially 16k 
(top left), 32k (top right), 64k (bottom left) and 128k (bottom 
right) stars. At the top of each panel the number of passed two- 
body relaxation times is shown (Eq. and crosses indicate 
the time of core collapse and of cluster dissolution. 



stars of the first generation show a more extended distribution up 
to 40 pc. 

Progressively the second generation stars spread out due to 
dynamical encounters so their radial distribution extends. The 
middle panels in Fig. [T] and |2] show that even after 5 Gyr of evo- 
lution the two populations still retain different distributions. At 
this time, which corresponds roughly to one initial relaxation 
time for the whole cluster, the extent of the second generation 
stars reaches 35 pc, but they are still more strongly concentrated 
to the central part. 

The bottom panels in Fig. [T] and |2] show that after nearly 
12 Gyr of evolution (slightly more than 2 initial relaxation times) 
the two populations have similar radial distributions and can no 
longer be distinguished on the basis of their dynamical proper- 
ties. 



3.2. Number of escaping stars 

As already stated in the introduction, two competitive processes 
act in the clusters: the loss of stars from the outer cluster parts 
will first reduce the number of bound first generation stars; and 
the dynamical spread of the initially more concentrated second 
generation stars will stop this differential loss when the two pop- 
ulations are dynamically mixed. We now describe the operation 
of this interplay. 

The time evolution of the half-mass radii of the two popu- 
lations is displayed in Fig. [3] for clusters with different initial 
numbers of particles. Stars which have already been lost by the 
clusters are not taken into account, so this radius reflects only 
bound stars. Initially, stars of the first generation have a larger 
half-mass radius as expected by our selection scheme (e.g. 5 
and 12 pc respectively for the case with = 128k stars). The 
fast initial increase of the half-mass radius is due to the evolu- 
tion of th e massive stars which lose a large part of their mass 
(Fig. 1 in lBaumgardt & Makinoll2003l) . Winds from these stars 
are not retained so the potential well of the cluster is lowered and 
the outer regions expand. Later the tidal forces of the host galaxy 
remove outer stars and the half-mass radius decreases until the 
final disruption of the cluster. The strong variation at the end of 
the evolution is due to the small number of low-mass stars still 
remaining in the cluster. They are also more pronounced in the 
smaller clusters. 

On the other hand, the second generation stars have initially 
a smaller half-mass radius as a result of their lower specific en- 
ergy. The initial expansion due to the evolution of massive stars 
has a smaller effect on the half-mass radius of the second genera- 
tion stars compared to first generation stars as second generation 
stars populate the centre and hence the deeper and less perturbed 
parts of the potential well. Later the half-mass radius increases 
continuously with time following the spread of the radial distri- 
bution as the cluster tries to achieve energy equipartition. After 
some time the two radii are similar and the further evolution of 
the two populations is identical. 

Fig.[3]also shows the number of passed two-body relaxation 
times, Nrei, defined as: 



A^rel(f) 



r' dt 
Jo ri-h(f)' 



(2) 



where t^hit) is the current relaxation time given by Eq. ([T) and f is 
the integration variable. The relaxation time, frh, decreases with 
time as both rh,n and N decrease (see Eq. ([TJ). For each clus- 
ter we also indicate the time at which core-collapse occurs and 
the time at which dissolution of the cluster happens {i.e. when 
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Fig. 4. Evolution with time of the ratio of first generation (dashed 
lines) and second generation (solid lines) low-mass stars to their 
initial content for the models with 8k (top left), 32k (top right), 
64k (bottom left) and 128k particles (bottom right). At the top 
of each panel the number of passed two-body relaxation times is 
shown (Eq. (|2]i), and crosses indicate the time of core collapse 
and of cluster dissolution. 
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Fig. 5. Number ratio between the second (with low initial spe- 
cific energy) and first (with high initial specific energy) popula- 
tion of low-mass stars in a cluster as a function of time for the 
same models as in Fig.|4] At the top of each panel the number of 
passed two-body relaxation times is shown (Eq. (|2]l), and crosses 
indicate the time of core collapse and of cluster dissolution. 



95% of the initial mass has been lost by the cluster) as com- 
puted bv lBaumgardt & Makinol (l2003b . These two times are la- 
belled by "cc" and "dis", respectively. Cluster homogenisation 
needs around 2 relaxation times nearly independently of the ini- 
tial number of stars. However there is still a weak dependence 
with the initial number of stars since the effect of mass loss 
does not scale wi th N in the same way as the relaxation time 
jBaumgardtll200 lb . Other radial properties, like the mean radius 
and the radial distribution of stars, are also similar in both pop- 
ulations when A^iei > 2. Since GCs are several relaxation times 
old, they have lost all information on their initial radial distribu- 
tion. The case of oj Cen, which is dynamically younger, will be 
discussed in § 3.4 

As previously seen, the effect of the external potential of the 
Galaxy on the cluster is to strip away stars with high specific en- 
ergy. Fig. ID shows the evolution of the fraction of stars in each 
population that still remains in the cluster as a function of time 
for different initial numbers of stars. Initially, only stars of the 
first generation populate the outer part of the cluster owing to 
their high specific energy. Therefore only first generation stars 
are lost in the beginning. This lasts until the second generation 
stars migrate into the outer part of the cluster Depending on the 
cluster mass, it takes between 1 to 4 Gyr to start a loss of second 
generation stars. This timescale scales roughly with the relax- 
ation time (see also § 3.3 for more details). Indeed, as shown 
before, more massive clusters have a more extended initial ra- 
dial distribution so the migration of the second generation stars 
lasts longer. Besides, these clusters also have a deeper potential 
well and the erosion by the tidal forces of the Galaxy is reduced. 

Due to the time-delay to lose second generation stars, their 
remaining fraction in the cluster is always higher than that of 



the first generation stars except during the final stage of cluster 
dissolution. Fig.|5]quantifies this point by showing the time evo- 
lution of the number ratio of second to first generation stars. As 
a direct consequence of our selection procedure, the initial ratio 
is about 0. 1 ; and it then increases gradually with time. It tends to 
stay nearly constant as soon as the two distributions are similar 
Finally, at the time of cluster dissolution, large variations occur 
due to the low number of low-mass stars present in the cluster 

Over the cluster history, the fraction of second generation 
stars relative to first generation ones increases by a factor of 2.5. 
Therefore, when the two populations have the same radial distri- 
bution, these stars represent 25% of the low-mass stars present 
in the clusters. Presently, observed ratios have been derived only 
spectroscopically for NGC 2808 (70%) and NGC 6752 (85%), 
since a large number of stars needs to be analysed. Thus, the in- 
ternal dynamical evolution and the dissolution due to the tidal 
forces of the host Galaxy are not efficient enough to produce the 
high observed value. An additional mechanism is thus needed to 
expel the first generation stars more effectively. 

3.3. Influence of the Initial number ratio 

The next point we want to investigate is how our results de- 
pend on the number of stars selected to be first or second gen- 
eration stars. Fig. |6] addresses this issue: for the computation 
with initially 128k stars, we select different number of second- 
generation stars which represent 1, 5, 10, and 25% of the initial 
mass of the cluster and follow how the ratio of second to first 
generation of stars evolves with time. For the 5% case, we re- 
trieve our previous results: the ratio increases by a factor of 2.5 
and then flattens until cluster dissolution. The more numerous 
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Fig. 6. Increase of the ratio of second to first generation stars Yig.7. Ratio of stars with negative to positive orbital angular 

as a function of time for different numbers of initially selected momentum in the cluster as a function of time. At the bottom, the 

stars for the model with N = 128k stars. 1, 5, 10, 25% refers cuiTently passed relaxation time is shown (Eq. ©), and crosses 

to fraction of the cluster mass populated with second generation indicate the time of core collapse and of cluster dissolution, 
stars. 



the second generation stars are, the stronger the increase is. In 
the case of 25% of stars selected, the ratio increases by a factor 
of 5. This trend can be understood as a result of the selection pro- 
cess which selects mainly those stars which remain bound to the 
cluster until the end. However this case can only be achieved if 
the IMF slope for the massive stars is around 1.7 (when Salpeter 
value is 2.35). 

For a standard IMF we expect only a limited number of low- 
mass stars to be second generation stars. Thus the models above 
with a low fraction of second generation stars are more realis- 
tic to what happens in real GCs and all of them predict only a 
modest increase of the fraction of the second generation stars. 

3.4. Angular-momentum-selected sub-populations 

In the previous sections we have only considered the case of 
two populations which initially diff'er by their spatial distribu- 
tion as we expect such diff'erences if a second generation of 
stars are born from the ejecta of massive stars. In this con- 
text, the GC cj Cen is an exceptional case as it not only shows 
variations in light elements but als o a wide spread in metal- 
licity from [Fe/H] = -2 to -0.5 (iNorris&Da Costal Il995h . 
Additionally, a close inspection of the color-magnitude diagram 
of a> Cen shows that several po pulations exist at all the evo- 
lutionary stages, from the RGB (Lee et al. 1999'; Pancino et al.' 
2000; Sollima et al. 2005a) to the sub-giant (Ferraro et al. 2004; 
SoUima et al.ll2005bt IVillanova et al. 2007j) branches and main 



sequences dBedin et al.ii2004: Piotto et aljl2005h . Interestingly, 



iNorris et all (119961) found the more metal-rich stars to be more 
spatially concentrated. Evidence for significant differences in 
the kinematics and the spatial distribution between the metal - 
rich and oth er stars have been collected (iNorris et al. 
|Jurcsikl[l998l) . However these findings have been challenged by 



iPancino et al.l ( l2007h on the basis of an analysis of a larger sam- 
ple. 

Therefore in this section we want to determine if the time 
needed for a cluster to lose its initial kinematic properties is the 
same as for the initial radial distribution. In particular we want to 
test how quickly the kinematic differences are erased since dif- 
ferences are observed in cj Cen whereas other GCs do not show 
such strong observational evidences. To assert this timescale for 
achieving kinematic homogenisation, we setup the initial con- 
ditions for the analysis in another way: we divide the low-mass 
stars according to their angular momentum about the cluster cen- 
tre along the z-axis (the cluster orbital motion occurs in the x-y 
plane with positive angular momentum). 

First, Fig.|2]shows the evolution of the number ratio between 
stars having negative and positive orbital angular momentum. 
Initially half of the low-mass stars populate the groups with 
prograde (L, > 0) and retrograde (L. < 0) motion. A rapid 
drop of the ratio of the number of stars with retrograde mo- 
tion (from 1 to 0.95) occurs due to the preferential loss of stars 
with negative L, during the evolution of the massive stars, fol- 
lowed by a very slow decline during the cluster evolution. At 
the time of core collapse, this ratio is around 0.9. This fast ini- 
tial change can be explained by the mass-loss of massive stars 
which lowers the potential well of the cluster leading to the ex- 
pulsion of the stars in the outer parts as we have seen before. 
This ejection is more pronounced for the stars with retrograde 
motion (L, < 0) relative to the cluster orbit than those with pro- 
grade motion {L, > 0). This finding seems to contradict pre- 
vious works which found that retrograde orbits are more sta- 
ble than prograde ones if the cluster is on a prograde orbit i n 
a steady external potential (see e.g. Fukus hige & Heggi3l2000l) . 
The difference can be exp lained by the different coordin ate sys- 
tems used in these studies: iBaumgardt & Makinol ( |2003|) applied 
an accelerated and non-rotating coordinate system (the Galaxy 
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(Gyr) 



Fig. 8. Fraction of stars which have kept their initial orbital an- 
gular momentum for negative (top) and positive (bottom) cases 
as a function of time. At the bottom, the currently passed relax- 
ation time is shown (Eq. (|2]l), and crosses indicate the time of 
core collapse and of cluster dissolution. 
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Fig. 9. (Top) Time needed for 1 % of the first generation (crosses) 
and second generation (full squares) of low-mass stars to be lost 
by the cluster as a function of the initial mass of the cluster. 
(Bottom) Ratio of these timescales to the initial relaxation time. 



rotates about the cluster) while an accelerated and ro t ating o ne 
is used by the computations of iFukushige & Heggid (|200C/, in 
the cluster centric frame the Galaxy is fixed; the cluster-centric 
frame is therefore non-inertial). R etrograde orbits in a rotating 
frame (iFukushige & Heggid 12000.) are found to be the stable 
ones, and a few of them shift to prograde orbits if we con- 
sider them in a non-rotating frame. This effect leads to the in- 
ve rsion of the dominant populat ion and explains the discrepancy 
of iBaumgardt & Makinol (l2003b with the older study. We con- 
clude that stars on retrograde orbits with respect to their cluster 
centre are preferably lost from the cluster, where prograde or- 
bits are in the same sense as the cluster orbit about the Galaxy. 
So an initially non-rotating cluster begins to rotate in a prograde 
manner after the massive stars have evolved. But the effect re- 
mains sm all (see Fig. |7]l. Similar results have been obtained by 
lYim & L ee (2002) who state that clusters with no initial rotation 
can gain a slight positive angular momentum (i.e. that the cluster 
rotates in the same direction that the orbital motion.). 

To measure the time needed to lose the initial angular mo- 
mentum distribution, we need to look more carefully at each 
population, as due to stellar encounters stars can change their 
orbital parameters and hence invert their orbital angular momen- 
tum. Initially only two populations exist: stars with L-{t = 0) < 
and L^(t = 0) > 0. As time proceeds, stellar interactions induce 
changes to the orbits of stars and some stars show an inversion 
in L,. Hence the cluster splits into four populations depending 
on their initial and current angular momentum. Fig.|8]shows the 
evolution of the fraction of stars which have kept their initial 
orbital angular momentum for the population of stars with nega- 
tive and positive L,. Here homogenisation will be obtained when 
50% of the stars have inverted their orbital angular momentum. 
Fig.|8]displays the fraction of stars keeping their initial AM sign 
rapidly decreases to 70-80%, thus 20 to 30% of stars suffer such 



an inversion during the early radial expansion following the evo- 
lution of the massive stars. Then the inversion rate decreases 
slowly and homogenisation fully occurs after around 2-3 relax- 
ation times. Thus the loss of information about orbital angular 
momentum follows the same timescale as the radial distribution. 

Therefore if there is any dynamical (radial or kinematic) 
property which differs initially between different cluster sub- 
populations, most clusters would have erased them during their 
evolution as they are several relaxation times old. cj Cen has the 
longest relaxation time which spans a range fro m several 10^ yr 
at the centre to 10'*' yr at the half-mass radius ( van de Ven et al.l 
2006). Thus it has evolved only for about one relaxation time and 
we should indeed see differences, at least at a low level, if the 
cluster was born with dynamically distinct populations. Thus if 
the results of Pancino et al. (2007) are confirmed we would have 
to conclude that the different populations observed in the CMD 
were born with the same radial and kinematic properties. 

3.5. Timescale for early processes 

Fig. |9] gives the timescales for a cluster to lose 1% of both gen- 
erations of low-mass stars as a function of the mass of the clus- 
ter This time corresponds to the end of the initial plateau seen 
in Fig.|4]for the second generation. Since IBaumgardt & Makinol 
( 1200 3) computed several N-body models with 8k and 16k stars, 
we have included all the data points. 

In these models, the first generation stars are the first to be 
noticeably lost after around 100 Myr and this is nearly indepen- 
dent of the cluster mass. On the other hand, the time until 1% of 
the second generation evaporates increases with the cluster mass 
because these stars have to travel to a larger radius before be- 
ing lost. For the heaviest clusters this happens only after 4 Gyr 
(this time scales as 0.7-0.8 times the relaxation time, see bottom 
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panel of Fig.|9l). It should be noted that the lightest clusters show 
a strong statistical effect due to the small number of second gen- 
eration stars initially present. As heavier clusters show an almost 
identical time relative to their relaxation timescale, we can quite 
confidently extrapolate our results to more massive clusters. 

As already explained, we need to find an additional mech- 
anism which would more vigorously eject first generation stars 
and at the same time keep the second generation stars. Several 
processes can play this role. First, disk or bulge shocking has 
been argued to heat up the outer part of the cluster when it passes 
through the Galactic plane. This effect increases the loss of stars. 
As the timescale between two passages in the Galactic plane is 
around 10** yr and the appearance of second generation stars in 
the outer regions needs a few 10^ yr, a high number of shocks 
affecting nearly exclusively the first generation should have hap- 
pened. Such process can also l ead to the destruction of many 
GCs over the Galactic history ( Shin et al] '2008). However old 
GCs formed about 12 Gyr ago whereas the Galactic disc formed 
continuously from 10 Gyr ago to the present. As we need an effi- 
cient physical mechanism which can act at the early time of GC 
evolution, disc shocking can be discarded to remove preferen- 
tially first generation stars. 

Initial gas expulsion by massive stars operates early in 
the cluster history (a Myr after cluster formation or earlier, 
iBau mgardt et al.' 2008b). As the residual gas still present after 
star formation is removed, the potential well of the cluster can 
be strongly reduced and the outer parts of the cluster can be- 
co me unbound. This p rocess has already been used successfully 
bv lMarks et al.l (l2008l) to explain the challenging correlation be- 
tween th e central conce ntration and the mass function of GCs as 
found bv lDe Marchiet al. (2007) . 

Gas expulsion is expected to proceed on a short timescale 
similar to the evolution-timescale of massive stars, and will take 
place well before any changes due to the internal dynamics or the 
tidal interaction with the host Galaxy, the two processes investi- 
gated in this paper However the results presented in the previous 
sections only apply for the limiting hypothetical case of a 100% 
star formation efficiency. 

Gas expulsion can occur at two different times: the intra- 
cluster gas can be washed out by the ionising radiation of OB 
stars, or later by the energy released by supernovae explosions. 
As these massive stars are also expected to eject the slowly es- 
caping matter from which the second generation stars form, we 
have two different chronologies to form a cluster with two chem- 
ically distinct populations of low-mass stars depending on the ef- 
ficiency of the UV radiation to eject the intra-cluster gas. Indeed, 
in the first case the slow winds of massive stars will accumu- 
late in the cluster centre where no pristine matter will be left 
while some dilution is expected to occur if the ionising radiation 
is inefficient to expel the gas from the cluster . This latter case 
would occur in clusters at the high-mass scale (iBaumgardt et al.l 
l2008bl) . In both cases the second generation stars have to be 
formed before the supernovae explode. 

However the ant icorrelation observed between Li and Na 
dPasquini et al.ll2005l) can help us to discriminate between these 
two possibilities. Indeed some Li is present in stars with abun- 
dance anomalies and this has to come from the intra-cluster gas 
itself as noted above. Therefore gas expulsion in GC that have 
abundance anomalies is more likely to be due to the supernovae 
than the UV radiation of m assive stars. Such models have been 
studied bv lGoodwinI d 1 997h . It should also be noted that the slow 
winds of massive stars will be released at the same time as the 
ionising flux from O stars, so that it would be difficult to expel 
the intra-cluster gas while leaving their slow winds untouched. 



We have thus found a physical reason as to why typically more- 
massive GCs show more pronounced abundance anomalies (see 
Fig. 12 in Carretta 2006) as these clusters are more able to keep 
their residual gas until r emoval though supernovae (see Fig. 3 in 
IBaumgardt et al.ll2008bl) . 

4. Conclusions 

We have analys ed the dynamical models of 

iBaumpardt & Makinol ( l2003h by selecting two stellar pop- 
ulations according to their specific energies with the aim to 
explain the high number of stars harbouring abundance anoma- 
lies in GCs. The basic idea that we followed is that stars with 
abundance anomalies stem from a second generation which 
formed centrally concentrated out of the slow winds of rapidly 
rotating O stars. We find the following results: 

- Stars with low specific energy are initially concentrated to- 
ward the cluster centre and progressively diffuse into the 
outer parts of the cluster as a result of the system attempt- 
ing to evolve to energy equipartition. About two relaxation 
times are needed to achieve a complete homogenisation be- 
tween the two populations corresponding to several Gyr of 
evolution. 

- At early times, stars with high specific energy are preferen- 
tially lost from the cluster due to the tidal forces of the host 
galaxy. Second generation stars are only lost after several 
Gyr when they have diffused into the outer cluster parts. 

- Any radial or dynamical differences between different stel- 
lar populations would be erased over the dynamical history 
of old GCs. Even with strongly different radial distributions, 
as we assumed in this paper, GCs would today present ho- 
mogeneous distributions. This finding could be in sharp con- 
trast with other explanations for the chemical enrichment of 
GCs. A scenario which relies on the en counters between 
red giant star s and main sequence ones d Suda et al.ll2007t 
lYamada et al.l [2008) should produce more second genera- 
tion stars in the cluster centre due to the high central den- 
sity. Unfortunately, we presently do not know the true radial 
distribution of stars with abundance anomalies in GCs. 

- The typical timescale for a cluster to lose information of its 
initial conditions on the angular momentum of stars is sim- 
ilar to the time needed for radial homogenisation. Thus all 
initial radial and kinematic properties of GCs are lost after 
about two relaxation times. 

- The number ratio between the second and first stellar genera- 
tion increases by a factor 2.5 durin g the cluster evolution (see 
also the results of lD'Ercole et al.l2008l) . The resulting ratio at 
old age would be too low to explain the ratio of anomalous 
to standard stars observed in GCs if we assume a standard 
IMF for the first generation stars. 

- The observed ratio would be achieved if the IMF of the first 
generation were top-heavy. This however would raise the 
problem of cluster sur vival (iBekki & Norrisll2006l) . 

- To reproduce the observed ratio for an invariant IMF, a 
more powerful mass-loss mechanism is needed. Residual 
gas expulsion due to stellar feedback from OB stars is 
the prime physical candidate occurring on a short enough 
time scale and mainly removing the outer part of the cluster 
(Mar ks et al.l 1200 8). Further investigations in this direction 
will be addressed in a forthcoming paper. 
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